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Abstract: In the context of Markov processes, both in discrete and continuous 
setting, we show a general relation between duality functions and symmetries of 
the generator. If the generator can be written in the form of a Hamiltonian of 
a quantum spin system, then the "hidden" symmetries are easily derived. We 
illustrate our approach in processes of symmetric exclusion type, in which the 
symmetry is of SU{2) type, as well as for the Kipnis-Marchioro-Presutti (KMP) 
model for which we unveil its SU{1, 1) symmetry. The KMP model is in turn 
an instantaneous thermalization limit of the energy process associated to a large 
family of models of interacting diffusions, which we call Brownian energy process 
(BEP) and which all possess the SU{1, 1) symmetry. We treat in details the 
case where the system is in contact with reservoirs and the dual process becomes 
absorbing. 
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1 Introduction 



Duality is a technique developed in the probabilistic literature that allows to obtain elegant 
and general solutions of some problems in interacting particle systems. One transforms the 
evaluation of a correlation function in the original model to a simpler quantity in the dual 
one. 

The basic idea of duality in interacting particle systems goes back to Spitzer [llj who 
introduced it for symmetric exclusion process (SEP) and independent random walkers to 
characterize the stationary distribution. Later, Ligget [8] systematically introduced duality 
for spin systems and used it, among others, for the complete characterization of ergodic 
properties of SEP, voter model, etc. Duality property might also be useful in the context of 
transport models and non-equilibrium statistical mechanics, that is when the bulk particle 
systems is in contact at its boundaries with reservoirs working at different values of their 
parameters. For instance, considering again the symmetric exclusion process in contact 
with particle reservoirs at different chemical potentials, Spohn used duality to compute 
the 2-point correlation function |12j . showing the existence of long-range correlations in 
non-equilibrium systems. In the case of energy transport, i.e. interacting particle systems 
with a continuous dynamical variable (the energy) connected at their boundaries to thermal 
reservoirs working at different temperatures, duality has been constructed for the Kipnis- 
Marchioro-Presutti (KMP) model [7] for heat conduction and also for other models [6]. 
Consequences of duality include the possibility to express the n-point energy correlation 
functions in terms of n (interacting) random walkers. Duality has also been used in the 
study of biological population models, see and references therein. 

One should notice that the construction of a dual process is usually performed with 
an ad-hoc procedure which requires the ansatz of a proper duality function on which the 
duality property can be established. The closure of n-point correlations functions at each 
order might be an indication that a dual process exists. However in the general case the 
closure property is neither sufficient nor necessary to construct the dual process. In this 
paper we present a general procedure to derive a duality function and a dual process from 
the symmetries of the original process. When applied to transport models, our theorems 
allow to identify the source of the existence of a dual process with the non-abelian symmetries 
of the evolution operator. The idea is simple: transport models have in the bulk a symmetry 
associated with a conserved quantity, the one that is transported. It may happen in some 
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cases that this symmetry is a subgroup of a larger group, i.e. that extra (less obvious) 
symmetry are present. In that case, one can describe the same physical situation as the 
transport of another quantity (another element of the group), and in some cases this makes 
the problem simpler. In the physics literature Sandow and Schutz [10] realized that this 
is case for the SEP process, whose SU{2) symmetry they made explicit by writing the 
evolution operator in quantum spin notation. In this paper we study in full generality the 
relation between duality and symmetries. We give a general scheme for constructing duality 
for continuous time Markov processes whose generator has a symmetry. For interacting 
particle systems used as transport models we detail the effect of the reservoirs. For particle 
transport models we generalize the symmetric exclusion process to a situation where each site 
can accommodate up to 2j particles, with j £ N/2. For energy transport models we uncover 
a hidden SU{1, 1) symmetry in a large class of models for energy transport (including KMP 
model) which explains their duality property, as the SU (2) does for the SEP process. 



2 Definitions and Results 
2.1 Generalities 

Let {rjt)t>o denote a Markov process on a state space Q. Elements of the state space are 
denoted by ^, C,-- The probability measure on path space starting from r] is called P,,, and 
denotes expectation with respect to P,,. In the whole of this paper, we will restrict to 
Feller processes. In that case, to the process {r]t)t>o there corresponds a strongly continuous, 
positivity-preserving, contraction semigroup At : C(J7) — > C(0) with domain the set C{0,) of 
continuous functions / : $7 ^ M 

Atfiv) := IE,/(r?t) = HfMlvo = v) = J f{v')Pt{v,dri') (1) 

where Pt{'r],drj') is the transition kernel of the process. The infinitesimal generator of the 
semigroup is denoted by L, 

t^o t 

and is defined on its natural domain, i.e. the set of functions / : — > M for which the limit 
in the r.h.s. exists in the uniform metric. We also consider the adjoint of the semigroup, 
with domain the set of signed finite Borel measures, : — > A4(0), defined by 

< > = < Atf,fi> 
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where the pairing < •, • >: C(0) x — > M is given by 

</>/">= y" fdl-t 

The processes which appear in our apphcations will always be either jump process or diffu- 
sions. 

Example 2.1. In the case that the Markov process {r]t)t>o is a pure jump process and the 
state space O is finite or countable then the generator is of the form 

Lf{r})=Y.c{r^,r^'){f{v')-f{v)) 

where c{i], rj') >0 is the rate for a transition from configuration ?y to configuration rj' . Equiv- 
alently we can write 

Lf{r,)=J2L{v,v')f{r]') 

where L{r],rj') is a matrix having positive off-diagonal elements and rows sum equal to zero, 
namely 

) c{r], T]') if rii^T]' 

In the context of a countable state space ft we have the usual exponential of a matrix, so that 



L{r],ri') 



At = e*^ = X;(tL)Vi! 



i=0 



and A1 = Af where the superscript denotes transposit 



ion. 

Example 2.2. General diffusion processes with state space O = are also considered 
here. In this case the generator take the form of a differential operator of the second order 

ij'=l ■' i=l 

(see 113^ for general conditions which guarantees that L satisfy the maximum principle and 
thus generate a positivity preserving semigroup). 
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2.2 Duality and Self-duality 

Definition 2.3 (self-duality). Consider two independent copies {r]t)t>o and {^t)t>o of a 
continuous time Markov processes on a state space We say that the process is self-dual 
with self-duality function Z) : O x Q — > M if for all (ry, ^) £ x Q, we have 

Er,D{r]t,0=^^Div,^t) (2) 

Definition 2.4 (duality). Consider two continuous time Markov processes: {r]t)t>o on a 
state space Q and {£,t)t>o on a state space ^duai- We say that {(,t)t>o is the dual of {rjt)t>o 
with duality function D : 0, x Qdual ^ if for all rj £ Q, £ ^duai we have 

If At denotes the semigroup of the original process {r]t)t>o and denotes the semi- 

group of the related dual process (^t)t>o then, using Eq. ([1]), the definition 12.41 is equivalent 
to 

AtD{r,,0 = Af^'''D{r,,() (4) 

where it is understood that on the l.h.s. of ([ID the operator At works on the rj variable, 
while on the r.h.s. the operator^d^""' works on the ^ variable. 

If the original process {i]t)t>o and the dual process {S,t)t>o are Markov processes with 
finite or countably infinite state space J7, resp. ^Iduai, (cfr. Example 12. ip property Q is 
equivalent with its "infinitesimal version" in terms of the generators 

In matrix notation, this reads 

LD = DL^^^i (6) 

where D is the matrix with elements D{ri,^) and {r],(,) G x ^Idual- Remark that in this 
case D is not necessarily a square matrix, because the state spaces and 0,fiual are not 
necessarily equal and or of equal cardinality. 

When $7 = flduai and At = Af^°'\ then an equivalent condition for self-duality (cfr. ([2|)) 

is 

LD = DL^ (7) 
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2.3 Duality and Symmetries 

We first discuss self-duality and then duality. We consider the simple context of finite 
or countably infinite state space Markov processes. In many cases of interacting particle 
systems, the generator is a sum of operators working only on a finite set of coordinates of 
the configuration. Therefore, showing (self)-duality reduces to showing (self)-duality for the 
individual terms appearing in this sum, which is a finite state space situation. 

Definition 2.5. Let A and B he two matrices having the same dimension. We say that A 
is a symmetry of B if A commutes with B, i.e. 

AB = BA (8) 

The first theorem shows that self-duality functions and symmetries are in one-to-one 
correspondence, provided L and are similar matrices, which is automatically the case in 
the finite state space context. 

Theorem 2.6. Let L he the generator of a finite or countahle state space Markov process. 
Let Q he a matrix such that 

L^ = QLQ^^ . (9) 

Then we have 

1. If S is a symmetry of the generator, then SQ^^ is a self-duality function. 

2. If D is a self-duality function, then DQ is a symmetry of the generator. 

3. If S is a symmetry of L^^ , then Q~^S is a self- duality function 

4- If D is a self-duality function, then QD commutes with L^. 

Proof. The proof is elementary. We show items 1 and 2 (item 3 and 4 are obtained in a 
similar manner). Combining ([9]) with ([8]), we find 

L{SQ~^) = {SQ~^)L^ (10) 

i.e., D = SQ~^ is a self-duality function (see Eq. d?])). Conversely, if D is a self-duality 
function, then combining ([9]) with d?]) one proves ^ for S = DQ. □ 
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Remark 2.7. Self- duality functions are not unique, i.e. there might exist several self-duality 
functions for a process. This is evident form the fact that if D is a duality function for 
self-duality, and S is a symmetry, then SD is also a duality function for self-duality. An 
interesting question is to study the vector space of self-duality functions, its dimension, etc. 
However this question is not addressed in this paper. See |^ for a discussion of this issue 
and some examples in the context of Markov processes with discrete state space. 

Remark 2.8. In the finite state space context, L and are always similar matrices \T^, 
i.e., there exists a conjugation matrix Q such that = QLQ~^ . In interacting particle 
system the matrix Q can usually he easily constructed. As an example, in the case that C 
has a reversible measure, i.e., a probability measure fx on Q such that 

fiiv)Liv,v') = ^<v')Liv',v) (11) 
for all rj, rj' G 0,, then a diagonal conjugation matrix Q is given by 

Q{v,v') = l^ir])6r,,n' (12) 
In general, if is a stationary measure then 

is the generator of the time-reversed process, which is clearly similar to I?" . Therefore, the 
similarity of L and is equivalent with the similarity of the generator and the time-reversed 
generator. 

Self-duality is a particular case of duality. To generalize Theorem 12.61 to the context of 
(general) duality we need the notion of conjugation between two matrices. 

Definition 2.9. Let A be a matrix of dimension m x m and let B be a matrix of dimension 
n X n. A and B are called conjugate if there exist matrices C of dimension rax n and C of 
dimension n x m such that 

AC = CB, CA = BC (13) 

We then have the following analogue of Theorem 12. 6[ 

Theorem 2.10. Let L and L^uai be generators of finite or countable state space Markov 
chains. Then we have the following. 
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1. If Q is the matrix that gives the similarity 

^dual = QLdualQ"^ (14) 

and C and C are the matrices giving the conjugacy between L and L^uai in the sense 
of definition \2.9[ then: 

(a) For any symmetry S of the generator L, D = SCQ~^ is a duality function. 

(b) If D is a duality function, then S = DQC is a symmetry of L. 

2. If Q is the matrix that gives the similarity 

= QLQ-^ (15) 

and C and C are the matrices giving the conjugacy between and L^ual sense 
of definition 12.91 then: 

(a) For any symmetry S of the transposed generator , Q^^SC is a duality function. 

(b) If D is a duality function, then QDC commutes with . 

Proof. The proof of item 1(a) is given by the following series of equalities 

L{SCQ~^) = SLCQ-^ = SCLduaiQ-^ = {SCQ~')L^,,i (16) 

The first equality uses the hypothesis of S being a symmetry of the generator L, the second 
comes from the conjugation of the generators, the third is obtained from the similarity 
transformation p4|) . If one recall ^ then Eq. (|16p shows that D = SCQ~^ is a duality 
function. The proof of the other items follow from a similar argument. □ 

3 Examples with two sites 

In this section we present a series of examples where particles jump on two lattice sites. We 
wish to show how (self)-duality can be established by making use of the previous theorems. 
To identify the symmetries we will rewrite the stochastic generator, or its adjoint, in terms 
of generators of some symmetry group. Some of the examples will be useful later for the 
study of transport models. In fact, many transport models such as the exclusion process 
have a generator that is written as the sum of operators working on two sites. 
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3.1 Self-duality for symmetric exclusion 



We first recover the classical self-duality for symmetric exclusion [8]. One has two sites 
(labeled 1,2) and configurations have at most one particle at each site. Particles hop at 
rate one from one site to another, and jumps leading to more than one particle at a site are 
suppressed. As usual we write 0, 1 for absence resp. presence of particle. The state space 
is then 0, = {00,01,10,11}. Elements in the state space are denoted as rj = (r?ir/2)- The 
matrix elements of the generator are given by Lqijo = -^^io,oi = 1 = — -Z^oi,oi = — -^^lO.iOi and 
all other elements are zero. 

To apply Theorem 1 2. 6 1 we need to identify a symmetry S of the generator. The transposed 
of the generator can be written as 



I2 



(17) 



where the operators J? with i G {1,2} and a G {+,—,0} act on a 2-dimensional Hilbert 
space, with basis |0) = (i), |1) = (?), as 



J' 



; 


: 




( 












V 


; 





(18) 



-1/2 \ 

and Ij is the identity matrix. The operators with a G {+,—,0} satisfy the SU{2) 
commutation relations: 



J? 



-2J" 



(19) 



from which we deduce (cfr (jl7p ) that L?^ commutes with the three generators of the SU{2) 
group, J'* = J" (8> I2 + li ® ^'^^ ^ ^ {+' ~' 0}- ^ possible choice for the symmetry of is 
then obtained by considering the creation operator J'^ and exponentiating in order to have 
a factorized form 



More explicitly, in the basis |0) ® |0), |0) (g) |1), |1) |0), |1 



7+ 1+ 
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|1), the matrix S is 
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1 
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We also need the similarity transformation between L and L . The matrix Q, relating L 
to its transposed, is the identity since L is symmetric. A duality function for self-duality is 
thus given hy D = Q~^S = S. Notice that D can also be written as 

^(??i^2,66)= n 

ie{l,2}:5i = l 

which is the usual self-duality function of [8]. 

3.2 Self-duality for 2j-symmetric exclusion 

Now we consider two sites with at most 2j particles on each site, with j € N/2. The state 
space is = ill X where Qi = {0, 1, . . . 2j}. The rates for transitions are the following: 
if there are i]i particles at site 1 and r/2 particles at site 2, a particle is moved from 1 to 2 at 
rate 771 (2j — 772) and from 2 to 1 at rate 7?2(2j — So in this case the generator is given by 

L{r]i7]2, r][r]2) = ?7i(2j - r?2)(5^i_i,,,^(5^2+i,^^ + mi^j - 7?i)(5^,+i,^'^(5^2_i,^^ 
-(7?i(2j - r/2) + ??2(2j - m))V,'?i'^'?2,»7^, 

The transposed of this generator can also be expressed as the scalar product between two 
spin operators satisfying the SU{2) algebra, namely 

= J+ (g, J- + Jf + 2 J? ® - 2fli (g) I2 (20) 

where the J", i G {1, 2} and a G {+, — , 0}, act on a (2j -|- l)-dimensional Hilbert space with 
orthonormal basis |0), . . . |2j) as 

Jth) = {2j-rj,)\rji + l) 

Ji\m) = ^ihi-i) 

JiW) = {m-3)\m) (21) 

The standard symmetric exclusion process of the previous section is recovered when j = 1/2. 
Reasoning as above, a symmetry of the generator is 

S = Si® 82 = e-^^ e^^ 

which has matrix elements 5'(ryir/2, ■^1^2) = 5*1(7/1,^1)52(7/2,^2) with 

S^iv^,C^) = {V^\e'^\C^)=C^~?') (22) 
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where we adopt the convention (^) = for m > n. 

To detect the matrix Q giving the similarity transform between L and (notice that 
L is not symmetric anymore for j 7^ 1/2) we make use of remark 12.81 and use the fact 
that the invariant measures of the 2j-symmetric exclusion process are products of binomials 
Bin{2j, p), with a free parameter < p < 1 (this will be proved in Theorem 14. 2p . Therefore, 
if we choose p = 1/2 then a possible choise is Q = Qi <^ Q2 with 

Q^ir]^,v'i) = ^,,,'^(^;|^ (23) 

Combining (|22p and (|23p . Theorem 12.61 then implies that a duality function for self-duality 
is given by 

D = Di(S)D2 = Q^^Si Q2 ^^2 

with 

D^{m,^i) = {Q7'Si){m,Cr) = (24) 

Later, in Theorem 14.21 we will give a probabilistic interpretation of this function. 
3.3 Self-duality for the dual-BEP 

This is a process that can be viewed as a "bosonic" analogue of the SEP (particles attract 
each other rather than repel with the exclusion hard core constraint). The state space is 
il. = ill X with rjj = N, i.e. we have two sites each of which can accommodate an 
unlimited number of particles. For iji particles at site 1, r]2 particles at site 2, the rate of 
putting a particle from 1 to 2 is given by 2r]i{2rj2 + 1) and the rate of moving a particle from 
2 to 1 is given by 2772(2771 + 1) . We will see later how this process arises naturally as a dual 
of the Brownian Energy Process (BEP), see Section [5] below. 
The matrix of the generator is given by 

Hmm,vW2) = 27/1(2772 + 1)^1 

+ 2772(27/1 + 1)(5,,/ ^^,^+1(5;,/ ^^2-1 
-(87/1772 + 2771 + 2r]2)6^^^ri[^ri2,ri'2- (25) 

The transposed of the generator can be written in terms of generators of a SU (1, 1) algebra 
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as follows. On each site i £ {1, 2} we consider operators with a G {+, —,0} given by 

K+lVi) = (7?^ + 1/2)1??. + !) 

K^\vi) = m\vi-'^) 

Kf\r^,) = (7?, + 1/4)1??,) (26) 
They satisfy the commutation relations of SU{1, 1): 

[Kf,K^] = ±Kt 

[Kr^Kf] = 2K^ (27) 
The transposed of the generator then reads 

= 4 (^K:^ ® + ® K+ - 2Kf ^K^ + ^li^l2^ (28) 

From the commutation relations, it is easy to see that L-^ commutes with K"^ = Kf ® I2 + 
li (8) K2, for a G {+, — , 0}. A possible symmetry is then given by the matrix 

S = 81(^82 = e^^ (S) 

which has matrix elements 5'(??ir?2, Ci'^2) = 'S'i(??ii ^i)'S'2(??2 5 ^2) with 

= <''■!'='' l«-> = pi^^ll|Ffp5^ (29) 

A similarity transformation = Q^^LQ to pass to the transposed is suggested (remark 
I2.8p by the knowledge of the stationary measure of the dual-BEP (see Theorem 15. ip 



Q(r?i??2,??'i??2) = Qi{m,Vi)Q2(,m,V2) 

with 



-1 



The self-duality function corresponding to 5 of (p9|) and Q of (|30|) then reads 

D{r]iV2,^i^2) = -Di(??i,^i)i:)2(??2,6) 

D^{V^,^^) = Q-\m,m)S^{m,^i) = 2^' (^_g.)|^(2g,-l)!! ^^^^ 
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3.4 Self-duality for independent random walkers 

This is a classical example which is included here for the sake of completeness. We have two 
site 1 and 2, and particles hop independently from 1 to 2 and from 2 to 1 at rate one. So the 
rate to put a particle from 1 to 2 in a configuration with rji particles at 1 and r]2 particles 
at 2 is simply r/i . The generator is given by the matrix 

A self-duality function is D = Di ® D2 with 

The invariant measures are product of Poisson distributions and a possible conjugation is 
thus given hy Q = Qi fSi Q2 with 

Qi{Vi,^i) = 

'h- 

As a consequence, a symmetry of the generator is given hy S = Si <^ S2 with 

= (AQOfa,^.) = . (32) 

This symmetry comes once more from an underlying structure of creation and annihilation 
operators satisfying the Heisenberg algebra. Indeed, if one defines for i £ {1,2} operators 
af and a~ which are represented on an Hilbert space with basis |0), |1), |2), . . . by operators 
working as 

o-ilVi) = VilVi-'^) (33) 
then one easily verifies the commutation relation 

K"'«i"] = ■ 

In terms of these matrices, the transposed of the generator reads 

L'^ = -{af I2 - li af){a^ I2 - li a^) (34) 

which commutes with a"*" = af I2 + ti af ■ The symmetry S in (j32p is then recognized 
as S* = ex.p{af) exp(a^). 
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3.5 Duality between independent random walkers and a deterministic sys- 
tem 



As an example of application of Theorem 12.101 we consider again a system of independent 
random walkers jumping between sites 1 and 2. We show that this system is dual to a 
deterministic system evolving according to ordinary differential equations. 
We consider the "abstract" operator L 

L = -{aX - a^){a{ - al^^) (35) 

where o^t operators satisfying the canonical commutation relations 

[a-,a+] = 4,l. (36) 

One way to represent the previous operator is by considering 

a,^ = a+ = Xi, iG{l,2} 

which obviously satisfy (j36p . In this case the operator (|35p takes the form 



dx\ dx2 

which is the generator of the deterministic system of differential equations 

dxi(t) 



dt 

dX2{t) 

dt 

with solutions 



-ixi{t) - X2{t)) 

{xi{t) - X2{t)) (37) 



_ Xi(0)+X2(0) ^ XijO) - X2{0) ^^2t 

2 2 

Another possible way to represent the operator ([35]) has just been seen in the previous 
paragraph. In this case the creation and annihilation operators are represented as matrices 
with elements given by (j33p and then the operator (|35p can be seen as the transposed of the 
generator for a system on independent random walkers 

Lduai = -(a| ® I2 - li ® a+) o {a{ ^ I2 - 1^ a^) . 
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It is immediately checked that the function 



D{x,C) = D{xi,Ci)Dix2,a) 

with 

D{x^,^^)=xf 
gives a conjugation between L and L^^^i-, namely 

LD{x,0 = DLl^M.O 
In this case, this relation reads more explicitely, 

= + nax^^+^x^^-^ - (ni + na)^^ 4^ (39) 

This implies that x{t) and ^(t) are each other dual with duality function D and the 
following relation holds 

a;i(t)«^X2(t)«2 =E5,,5,(a;i(0)«iWx2(0)«2W) (40) 

where the expectation in the rhs is over the independent random walkers starting from initial 
configuration with r/i particles at 1 and rj2 particles at 2. We will come back to this example 
in section 16.41 



Remark 3.1. In the last example, we can still use other representations of the operators 
a~,af, satisfying the commutation relation [a^,a^] = dij, such that the abstract operator 
psp is the generator of a Markovian diffusion process. E.g., if we choose 



d 

■T, 

OXi 



_d_ 

dxi 



(41) 



then the abstract operator (j35p reads 



d d , f d d 



^ ' dx\ 8x2) ~^ ^'^^ \dxi dx2 
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which is the generator of the (degenerate) diffusion: the "coordinate {xi —X2)/2 undergoes 
a Brownian motion and (xi +X2)/2 remains constant. So in that case we also have duality 

^x^,x2{D{xi{t),ni)D{x2{t),n2)) = E„,^,„2(L>(a;i, ni(t))D(x2, n2(t))) 

where D{x,n) can he found by the recursion 

D{x,n + 1) = a'^D{x,n) 

e.g. the first five polynomials are 

D{x, 0) = 1, D{x, 1) = X, D{x, 2) = - 1, D{x, 3) = -3x + x^ D{x, 4) = 3 - + x'^ 

4 Symmetric exclusion processes 

In this section we study the 2j-SEP (with j S N/2), i.e. exclusion processes with at most 
2j particles per site, on a graph S. We show how we can understand self-duality for the 
2j-SEP from "classical duality" (in the sense of [8]) of the symmetric exclusion process on 
special graphs. We also consider two limits j — > oo leading to a deterministic process or a 
system of independent random walkers. Finally, we consider the boundary driven case, and 
show that we have a dual with absorbing boundaries. 

4.1 Symmetric exclusion on ladder graphs 

Consider a countable set 5, to be thought of as the underlying lattice, and a finite set 2 with 
cardinality 2j G N. The set I is to be thought of as a "ladder" on each site with 2j levels. 
The state space of SEP on the ladder graph S x I is Q, = {0, l}*^^-^. A configuration 
S is called finite if it contains a finite number of particles, i.e., if X]je5 aei C(^' ^) ^ 
The process is described as follows. Let p{i,l) denote a symmetric random walk kernel on 
S, i.e., p{i,l) = p{l,i) > 0, Z^/g5P(«,0 = 1- each site i £ S and level a £ I, there is at 
most one particle. Each particle attempts to jump at rate p{i, I) to a site / G S and level 

pel. 

More formally, the SEP on a ladder graph 5 xX is the process with the following generator 
on local functions / : O — > M 

Lf{c) = Yl E pii^mhc.){i-ai,mfi&''^'^''^^)-fic)) m 
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where ^(*'°)'('''^) denotes the configuration obtained from r] by removing a particle at site i 
level a and placing it at site / level (3. Since this process is a symmetric exclusion process 
on a special graph, then it is self-dual in the following sense: 

Proposition 4.1. Define for C,, C, & Q., C, finite, 

then we have the duality relation from 

E^DiCtX)=^^D{CXt) (43) 



where Ct, Ct cLf^ independent copies of the ladder SEP with generator (j42p starting from C,, 
resp. C,. 

4.2 From the ladder SEP to the 2j-SEP 

To define the 2j-SEP on a graph S we consider, for a given SEP on a ladder graph S x I 
with 2j levels, the process which consists of giving at time t > 0, and each site i G 5 the 
number of levels {i,a) which are occupied in Q. More precisely, define the map vr : ^ 
= {0,l,...,2j}'5 

C ^ vr(C) = V with T]i = Y, a) (44) 
Then we have the following theorem. 

Theorem 4.2. Let Q be the ladder sep with generator ()42p . Then the following holds: 
0-) r]t = '^iCt) is a Markov process on Q^^^ with generator 

L^'^fiv) = P(^'0^..(2i - m)if{r]''') - f{r])) (45) 

This process will be called the 2j-SEP or reduced ladder SEP with 2j levels, 
h) The process rjt with generator L^^^ is self-dual with duality function 
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for < rj a configuration with a finite number of particles (D is defined to he zero in 
other cases). More precisely, we have 

E.,D{rit,0=^iD{r,,Ct) (47) 

c ) The extremal invariant measures are of the form 

=®i65Bin(2j>i) 

where pi is harmonic for p{i,l), i.e., such that 

^p{i,l)pl = pi 
I 

In particular, if the only hounded harmonic functions are constants, then the only 
extremal invariant measures are products of hinomials with constant density. 

Proof. For point (a) remark that the jump rates in the generator (|45p only depend on the 
number of particles at a site, and not on the levels. Therefore, if / : — > M depends on C, 
only through rj = vr(C), i.e., if /((") = ^p{7^{C)) = ip{r]), then 

L/(C) = L(^V(77) (48) 

Therefore, for every local function : 0'--'^ — > R, 

V'C^CCt)) - ^(vr(Co)) - f L'^^\^){7:{Qs))ds = Mt (49) 

JO 

is a martingale w.r.t. the filtration Tt = <^{t^{C)s : < s < i}. This shows that ■K{Ct) is 
a solution of the martingale problem associated to the generator L^^\ and hence coincides 
with the unique Markov process generated by L*^-') (see Th. 4.1, page 182, of [1]). 

Now we turn to point (b). At each site i £ S we choose levels at random. For a 
given configuration r] G Qj, we choose £ Q consistent, i.e., such that 7r(^) = r]. Then the 
probability (w.r.t. to the random choices) that all chosen levels are occupied in ^ is exactly 
equal to D[r],^). As ■^{Ct) = ^t; the probability that the chosen levels are occupied at time 
t (i.e., in C,t) is given by ¥.jjD{r]t,^). By self-duality of the ladder SEP (Prop. 14. Ij) . the event 
that at time i > the chosen levels are occupied is the same as the probability that the 
particles evolving from the chosen levels during a time t find themselves at positions which 
are occupied in The latter probability equals E,^D{^t,il)- 
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Point c) follows from the fact that for the ladder SEP with generator (j42p . the product 
Bernoulli measures indexed by harmonic functions of are the extremal invariant mea- 

sures ( see [8] for details) and the image measure of a product of Bernoulli measure under tt 
is a product of Binomial measures. □ 

4.3 Limiting processes as j — > oo 

In this section we show that for large j the 2j-SEP converges, when considered on an ap- 
propriate time scale, either to a system of independent random walkers or to a deterministic 
limit, depending on the initial density. We remind the reader that for independent random 
walkers on a graph S, the configuration space is Ooo = N*^ and the generator is given by 

L^r^f{v) = E V{i)p{i,l){f{v''') - fm (50) 

The stationary measures are products of Poisson measures, and the process with generator 
(150p is self-dual with self-duality function 

A™(,.«=n^ (51) 

for finite configurations ^ < ?/, and D = otherwise. 

The relation with the reduced ladder SEP for large j is given in the following theorem. 

Theorem 4.3. Consider the process {rjl''^ ■ t > 0} with generator (I45|) started from initial 
configuration rj^^^ G Suppose that, as j oo, j/-'^ rj G ^oo, then the process 

{Vtj^j • ^ — 0} converges weakly in path space to a system of independent random walkers 
with generator (j50p and initial configuration rj. 

Proof. The process {fllj^j • ^ — 0} generator 

4 = ^L^'^ (52) 

In order to have a sequence of processes all defined on the same sample space we 
consider the auxiliary process on N"^ with generator 

^'jfiv) = E P('^ (l - 1^) ^('^^ ^ 2i)/(7?; < 2j) [fit') - /(r?)) (53) 
i,i&s ^ -^^ 
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This auxiliary process behaves as the process with generator Clj except for sites which have 
more than 2j particles, which are frozen. Started from an initial configuration with all sites 
having at most 2j particles, this process coincides with the process ??t/2j- For any local 
function / : N'^ M, we then have 



hm d'^fiji) = hm Pih l)m (l - |:) I{m < 2j)I{m < 2j) (/(r?*'') - /(t?)^ 

= Cirwfiv) (54) 



Therefore, by the Trotter-Kurtz theorem (see Theorem 2.5 of [4]), this implies that the 
corresponding processes ??t/2j converge weakly on path space as j — > co to the process with 
generator Cirw D 

To see that (jSip is (up to a multiplicative consant) a limit of duality functions of the 
2j-SEP, we start from the symmetry (j22p and use the similarity transformation with Q^^^ 
given by 

with 



Then the duality functions 
with S defined in (j32p . satisfy 

lim D^^^ = eDirw ■ 

j->oo 

Another possible limit is obtained when the initial condition has a number of particles 
that diverges proportional to j. This limit, as can be understood from the law of large 
numbers, is deterministic. 



Theorem 4.4. Consider the process 
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and suppose that (0) Xi{0) G [0, 1] for all i ^ S as j ^ oo. Then we have that x\ (t) 
converges to a deterministic system of coupled differential equations with generator 




Proof. The generator of the process x^ {t) reads 

i,i 

where x^-'-''*'^ arises from x^^^ by removing a unit l/2j from i G 5 and putting it at / G S. 
Therefore, for a local smooth function / : [0, 1]"^ ^ M we have, by Taylor expansion 




The result then follows once more from the Trotter Kurtz theorem. Since the limiting 
generator is a first order differential operator, the corresponding process is deterministic. □ 

4.4 Boundary driven case 

We first discuss a duality theorem for standard (i.e. j = 1/2) symmetric exclusion with extra 
creation and annihilation of particles at the boundary. The context is a countable set S, of 
which we distinguish a subset dS C S called the boundary. We then consider the generator 

^fiv) = E^'(^'0(^«(i-^(0)(/(^''')-/(^)) 

+ (1 - p,Mi) {f{if) - fiv)) + P.(l - Vii)) (/(V) - fiv)) (56) 
ieds 

where < pi < 1 represent the densities of the particle reservoirs with which the system 
is in contact at the boundary sites, and where ry* is the configuration obtained from r] by 
flipping the occupancy at i. 

The first part of the generator represents the hopping of particles on S according to a 
symmetric exclusion process, whereas the second part represents creation and annihilation 
of particles at the boundary sites. 

To introduce duality for this process, we introduce a set deS of sink sites, and a bijection 
i if. which associates each site i G dS to a sink site. The dual process will then be a 
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process that behaves as the exclusion process in the bulk, but particles can leave the system 
via boundary sites to sink sites, and will then be stuck at sink sites. More precisely, a 
configuration of the dual process is a map 

such that (^{i) G {0, 1} for i G 5 (only the sink sites can contain more than one particle). 
We call i^duai the set of all configurations of the dual process. 

The duality function is defined as follows: for r] G {0, l}"^, ^ G ^duai 

D{v,o= I n I n (^7) 

Ve5:5(i)=l / \ied,S J 

The idea here is that we have the ordinary duality function for the sites i £ S and for the 
sink sites, we replace the variable r]i by its expectation pi, corresponding to the stationary 
measure of the boundary generator Lj. 

The generator of the dual process is then defined as follows: 

Cduaifio = E p^'^ - m){fie'') - /(e)) + Yl ^(^) (/(^''^) - /(^)) (^s) 

i,l£S iedS 

We then have 

Theorem 4.5. The boundary driven exclusion process {'i]t)t>o with generator C in (|56p and 
the process {S,t)t>o with generator Cduai 'in (j58j) are dual with duality function D{r],(^) given 
by dnZD, i.e., 

Er,D{iit,0=^iD{ii,^t) (59) 

Proof. Abbreviate 

Lif{7j) = (1 - p,)rjii) {f{if) - firj)) + p,(l - 7/(0) (/(V) - fiv)) (60) 

and 

Lf'''f(.o = m{f{e^'')-f{0) (61) 

For f{r]) = r]{i) one sees that 

Lifiv) =Pi- Vii) (62) 
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and hence for ^ a dual configuration which is non-zero only on the sites i £ dS and on the 
corresponding sink site ie G deS, we find 

= mp'r{p^-vii)) 

From that and the self-duality of the symmetric exclusion process, it follows that 

CD{i^,0 = CduaiD{v,0 (63) 

□ 



In order to apply this duality result for the boundary driven process with generator (psj) , 
we first look at the boundary driven exclusion process on a ladder graph. More precisely, 
for ( G {0, we consider the process {Ct)t>o with generator 

^/(C) = E E P(^'0(C(^,«)(l-C(/,/3))(/(C(*'")'(''^^)-/(C)^ (64) 

+ E E(i - ^^)c(^' «) (/(^^''"^) - f(o) + c{i, «)) (/(c^^'")) - /(c)) 

iedS ael 

In words, this process is the ladder SEP, with additional boundary driving, where the creation 
and annihilation rate of particles at the boundary sites does not depend on the level. If we 
consider the reduced process, consisting of counting at each site i G S how many levels in 
I are occupied, i.e. r]t = T^iCt) then we recover once again a Markov process (cfr. Theorem 
14. 2p . This process, defined on the state space Q^^^ = {0, 1, . . . 2j}^ and called the boundary 
driven 2j-SEP, has generator 

Cjfiv) = E?'(^'^)(^(^)(2i-'?(0)(/(^*'')-/(^)) 

i,l€S 

+ J] (1 - p,)v{i) (/(?7(*'")) - fir])) + p,,(2j - rj,) {f{r^') - /(r?)) (65) 
ieds 

It turns out that the boundary driven 2j-SEP has a nice dual too. To introduce this dual, 
we consider admissible dual configurations as maps ^ : 5 U dgS — > N such that < ^(i) <2j 
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for i G S (only the sink sites can contain more that 2j particles) . The generator of the dual 
of the boundary driven 2j-SEP is a process on admissible dual configurations, defined by 



cf^'fio = ^Pi^,mmj-mmie'')-fiO) 

Prom Theorem 14.51 we then infer, in the same way as we derived Theorem 14.21 the following. 

Theorem 4.6. Let {rjt)t>o denote the boundary driven 2j-SEP with generator ()65p . Then 
{rjt)t>o is dual to the process {£,t)t>o with generator ()66p with duality function given by 

D{V,0= UpfUw^ (67) 

Proof. Denote for i £ dS, 

Cifiv) = (1 - p^Hi) (/(r?(''")) - fiv)) + P.(2j - Vii)) (/(V) " fiv)) 

and 

/:f'''/(0=e« (/(r^)-/(e)) 

One then easily computes that for a dual configuration which is non-zero only on the sites 
i £ dS and on the corresponding sink site ig S deS, 



(Vi 



The result then follows from combination of this fact and the duality relation (|47p . □ 



5 The Brownian Momentum Process and SU{1,1) symmetry 

In this section we study the Brownian momentum process that was introduced in [5l [2]. 
We will recover duality [6] in the context of our main Theorems and study the reversible 
measures of the dual process. 
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5.1 Generator and Quantum spin chain 

Let 5 be a countable set and p{i,j) a symmetric random walk transition probability on S. 
The Brownian momentum process is a Markov process {xt)t>o on M.'^ , with generator 

L=Y,p{i,j)Uj (68) 

where 

f d d V 

^'^• = V"^a^"^^'9^J ^^^^ 

and p{i,j) is a symmetric random walk kernel on S. 

The generator Lij conserves the energy xf + x'j and generates a Brownian rotation of 
the angle 9ij = aictan^Xj/xi). The interpretation of the generator (j68p is then as follows: 
each bond independently, at rate p{i,j) undergoes a Brownian rotation of its angle 9ij = 
arctan(xj/xi). An important example to keep in mind is 5 = Z"', and p{i,j) the nearest 
neighbor symmetric random walk. 

The processes with generator L can be related to quantum spin chains [6]. Consider the 
operators 



1 



2 



2 

1^ 

2 5x? 



which satisfy the commutation relations of SU{\, 1): 

[KlKf] = ±Kf 

[K7,Kt\ = 2Kf (71) 

Then the negative of the adjoint of the generator L can be seen as the quantum "Hamilto- 
nian" 

H = -L* = -4Y^ p{i,j) (k+K- + K-K+ - 2KfK^ + i") (72) 
with spin satisfying the 5C/(1, 1) algebra (in a representation with spin value 1/4). 
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5.2 Dual process 

In [6] we showed that the process with generator L in (j68p and (j69p has a dual, which is a 
system of interacting random walkers on S. We show here how this dual process comes out 
of the structure of the Hamiltonian (|72p . 

We notice that the SU (1,1) group admits a discrete (infinite dimensional) representation 
as (unbounded) operators on /2(N): 

/c+ie.) = Q + e.)lc^ + l) 

= + (73) 

where i S 5 and .^j £ N and |0), |1), |2), . . . denotes the canonical basis on l2{^)- It is imme- 
diately checked that the (unbounded) operators in (j73|) satisfy the 5C/(1,1) commutation 
relations in (j7ip . We then define a new generator via the same Hamiltonian as in (j72p . but 
now in the representation ([73]): 

From the previous equation and using the representation (i73]l we deduce that the Hamilto- 
nian above defines a Markov process {£,t)t>o with state space N*^ and generator 

where 

Lfr'fio = 2^.(2c, + i)(/(r'j)-/(e)) + 

2^j{2^i + l){f{en-f{0) (76) 

with ^^'-^ the configuration obtained from by removing a particle from i and putting it at j. 
Note that, in general, changing a representation does not imply that a generator continues 
to be a generator: the fact that H and Hfi^ai are well-defined as a Hamiltonian is conserved 
by similarity transformations (change of representation) but their property of being (minus) 
the adjoint of the generator of a Markov process is dependent on the representation, and 
needs to be verified by hand. 
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5.3 The duality function explained 

In [6] we found that L and L^uai are dual processes, with duahty function 

D{x,i) = J\Di{x,,i,) (77) 

with 

A(x,,,eO=(^f2_ (78) 

We now show how these functions arise from the change of representation. Suppose that we 
find a function 

C = C{x,i) = J{C,{xi,i,) 

such that 

K^Ci = CilCf (79) 

for a G {+, — ,0}, i £ S and Kf-, resp. /C? defined in ([70]) . ([73]) . The "matrix product" in the 
Ihs of (j79p is defined as the differential operator working on the Xj- variable of Ci{xi,(,i), 
and in the rhs amxi,^^) = C(xi, eD/C^lCUi)- 

Then for the generators L in ([68]) . ([69|) and the generator Lfi^ai in (f75|) . ([76]) we find, as 
a consequence of (|79p and using that L* = L, 

LC = L*C = CLl^^i (80) 

i.e., such a function C is a duality function (cfr. 

The equation (|75]) is most easy of a = +, it then reads 

^x?a(xi,ei) = + Ci(x„e, + 1) (81) 

To find Ci{xi,Q) we use KfCi{xi,0) = (that follows from ([73]) . ([79]) ) so we can choose 
Ci{xi,0) = 1 and then find, via (|81|) 

a(x,.{,) = ^^^ (82) 



which is exactly the duality function that we found in [6]. It is then easy to verify that ([79 
is also satisfied for a G {— , 0} with the choice 
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5.4 Reversible measures of the dual-BEP 



The dual of the BEP, with generator Lduai in (|75p and (j76p . is in itself an interacting 
particle system (particles attract each other), and it can therefore be considered as a model 
of independent interest. In some sense, it can be viewed as "the bosonic counterpart" of the 
SEP. Surprisingly, despite the interaction, the process has reversible product measures, as 
is shown below. Remark that, due to the attractive interaction between the particles, this 
process does not fall in the class of "misantrope processes" considered in [1], [3] (where one 
also has in particular cases stationary product measures, despite interaction). 

Theorem 5.1. Consider, for A < 1/2 the translation invariant product measure vx on 
with marginals 

Mm = k) = ^^^^^A'^ (83) 

where the normalization is 

Then vx is reversible for the process with generator L^^ai (|75p and (j76p . 

Proof. From the generator ([75]) . ([75|) . we infer that 

aiaj2i{2j + 1) = 2(j + l)(2i - l)aj+iai-i 

is a sufficient condition for detailed balance of a product measure with marginals /i(r/o = 
k) = for the generator Lf'^"'\ which is sufficient for detailed balance for L^uai- This leads 
to 



Oj ~ 2j + 2 
for some positive constant c. This in turn gives 



(2c) 



which is ()83p . The explicit expression for Zx can be obtained easily from the identity 

/oo 
^2fcg-xV2^^^^(2A;-l)!! 
-oo 

□ 
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6 The Brownian Energy Process 



As it was done for SEP, it is interesting to consider the Brownian Momentum Process on 
ladder graph S x I with |T| = m G N levels and look at the induced process which gives the 
energy at each site. 

6.1 Generator 

Consider the generator, working on local functions / : M"^^-^ — > M 

m. 

L=Y1 ^ Pih3)Lia,j(3 (85) 
jjG5Q,/3=l 



d d 



J2<^ (87) 



where 

In this section we show that for the process with the generator above, the total energy per 
site defined via 

m 

Zi = 

a=l 

is again a Markov process 

Theorem 6.1. Consider the process x{t) = {xi^a{t))i^s,a=i,...,m with generator L of (|85p 
and ()86p . Consider the corresponding process z{t) = {zi{t))i^s defined via the mapping (1870 . 
This is a Markov process on with generator 

L(-) = E^'(^'J-)4f (88) 

ij&S 

with 

and with stationary measures which are product measures with chi-squared marginals. 

Proof. Denote vr : {xi^a)i^s^a=i,...,m ^ {zi)i^s- Denote by di partial derivative w.r.t. zi and 
by di^a partial derivative w.r.t. xi^a. Then, using the identities 

= 2d, + 4xj^df 
di,adj,i3 = 4:Xi^aXj^pdidj (90) 
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we find for a function / : M"^^-^ — > M depending on x only through z = 7r(x) 

L/o^(x) = (L(-)/)(7r(x)) (91) 

The proof then proceeds via the martingale problem as in the proof of Theorem 14. 2[ The 
stationary measure of the process with generator L^'*") is deduced from the knowledge of the 
stationary measure for the process with generator L. Indeed, it is easy to check that for the 
process Xi^a{t), products of Gaussian measures ®iG5i,a^(0, cr^) are invariant and ergodic. 
The image measure under the transformation 7r(x) = z are products 'S'ieSiXmi^) where for 
£7 = 1, Xm(l) is the chi-squared distribution with m degrees of freedom, and for general a, 
follows from scaling Xmi^'^) = <^^Xm(l)- ^ 



6.2 Duality 

In this section we show that the BEP defined above has a dual process which is again a jump 
process. 

To construct the dual we follow a procedure similar to the one of the previous section. 
Remark that the generator L in (j85p and (j86p can be written in terms of the operators 



m 

K<-) = Y^Kl^ (92) 



with a G {+, — , 0}, where 



1 a2 



K, 



2 5x? 



In other words L is related to a quantum spin chain with Hamiltonian 

(94) 
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where the i^'-operators in (j92p and (j93p satisfy the commutation relations of SU{1, 1). More- 
over the K operators defined in ()92p . (1930 . can be rewritten in z-variables as follows: 



I 2 * 

i^o.M = ,^d^ + ^ (95) 

The generator L^™) of (|88p is then simply minus the adjoint of the Hamiltonian in 
(|94p. rewritten with i^-operators in z- variables. 

At this point it is important to remark that the SU{\^ 1) group admits a family of discrete 
infinite dimensional discrete representation labeled by m G N given by 



1 



K^T '\^^) = (j+^^jie.) (96) 

where i € S and S N and |0), |1), |2), . . . denotes the canonical basis on hiN)- We then 
define a new generator via the same Hamiltonian as in (j94p . but now in the representation 



(97) 

Using the representation (|96p we deduce that the Hamiltonian above defines a Markov process 
{it)t>o with state space N*^ and generator 

where 

^M,d.a/^(^^ = 2e.(2e, + m)(/(r'^) - /(o) + 

2^,(2^, + m)(/(e^)-/(0) (99) 
The duality function are given in the following theorem. 
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Theorem 6.2. The processes with generator L^"*) and L^^^i are each others dual, with 
duality function 

D{z,i) = J{Di{z,,^i) (100) 

where 

D^{z^M)=zf-—}i^ (101) 



with r(t) = /q°° a;* ^dx the gamma function. 



Proof. Let 

Ci{zi, ^i) = z, 



r (f ) 



One verifies easily that 

for a = +, — , 0. The proof then continues as in section [531 □ 

6.3 The instantaneous thermalization limit and the KMP process 

In the KMP model, introduced in [7], the energies Ei of different sites i G S are updated 
by selecting a pair of lattice sites (i,j) and uniformly redistributing the energy under the 
constraint of conserving Ei + Ej. In this section we show that the KMP model arises by 
taking what we call here an instantaneous thermalization limit of the process with generator 
L^^\ for the case m = 2. 

We start by computing the stationary measure of the process with generator L^^^ . 

Lemma 6.3. Let [zi(t), Zj[t)) be the Markov process with generator L^ij^\ starting from 
an initial condition (zi(0), Zj(0)) with Zi{{)) + Zj(0) = E. Then in the limit t — > cxo the 
distribution of {zi{t), Zj{t)) converges to the distribution of the couple {{E + e)/2, {E — e)/2) 
where e has probability density 

/(6) = C™(i?2_e2)f-i (102) 
—E < e < E and / = otherwise, and where Cm is the normalizing constant. 
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Proof. Define {E{t), e{t)) = {zi{t) + Zj{t), Zi{t) — Zj{t)). Then simple rewriting of Ll™' in the 
new variables yields that (£'(t),e(t)) is a Markov process with generator 

C' = 4(^2 _ ^2^q2 _ 4^gg^ (;^Q3^ 

From the form of C' we see immediately that E is conserved and that for given E, e{t) is an 
ergodic diffusion process with stationary measure solving 

d^{4{E^ - €^)f) + 5e(4me/) = (104) 

Now notice that the rhs of p02|) solves 

a,(4(i?2 - e^)f) + [Amef) = 

and hence p04|) . □ 

Denote by 7^ the distribution of {{E + e)/2, {E — e)/2). We can now define what me 
mean by instantaneous thermalization. 

Definition 6.4. Let f : [0,00)*^ M. For e = {ei)i^s a configuration of energies, {i,j) S 
S xS, {e'-,e'j) G [0,00) X [0, 00) we denote by t{e, e[, e'^) the configuration obtained from e by 
replacing Cj by e[ and ej by e'j. The instantaneous thermalization of a pair {i,j) £ S x S is 
defined by 

^rV(e) = I f{t{e,ele'^))djm{e'„e'j) (105) 

The instantaneously thermalized version of the Brownian energy process is then defined 
as the process with generator 

^ilfie) = PihJW^^fie) - /(e)) (106) 

This means that with rate p{ij) a pair {i,j) £ S x S is chosen and the energy is instanta- 
neously thermalized according to the measure 7^. From (jl02p one sees that, for m = 2, the 
uniform redistribution of the KMP model is recovered. 

It is interesting to consider the dual of the instantaneous thermalization process for 
general m G N. From the previous discussion one knows that in the case m = 2 this is 
just the dual of the KMP model. However, the model with generator has for general m a 
dual with different duality functions as is shown in Theorem 16.61 below. To introduce this 
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dual, we remind the reader that the Brownian energy process with generator L^"^^ is dual 
to the discrete particle jump process with generator L^J^^i- The following lemma gives the 
stationary measure of the dual BEP, which is needed in the construction of the instantaneous 
thermalized version of the dual BEP. 

Lemma 6.5. Let {kt,lt) evolve according to the generator C^J^^i, and suppose that initially 
ko + Iq = N, then in the limit t — > oo, (kf, It) converges in distribution to {{N + A)/2, (A^ — 
A) /2) where A has distribution ^ on {—N, —N + 2, . . . , A^} with 

/i(A) _ (iV-A + l)(jV + A-l + m) 
/x(A-2) {N + A){N-A + m) ^ ' 

In particular for m = 2, (ktJt) converges to the uniform measure on the set {{a,b) S 
{0,...,N} :a + b = N}. 

Proof. The process {Nt, At) := {kt + It, kt — It) performs transitions (N, A) (N, A — 2) at 
rate j{N + A){N - A + m) and (iV,A) (iV, A + 2) at rate \{N - A){N + A + m). The 
marginal At is then an irreducible continuous-time Markov chain on the set {— A^, —N + 
2, . . . , N}, and hence has a unique stationary measure. Since it is a pure birth and death 
chain, this measure is also reversible. The recursion (|107p then follows from detailed balance. 

□ 

We denote by "ym{k,l) the stationary distribution of lemma (j6.5p . For ^ G N"^, and 
{i,j) G 5 X 5, iCi^Cj) G N X N we denote by the configuration obtained from ^ by 

replacing the value at i by ^- and at j by . We then define the dual thermalization by 

^r''^"^V(0= E fm,e^,e,))lmie^,e,) (108) 

and the dual instantaneously thermalized energy process as the process with generator 

^I'lr'^m = E c^r'^^^/io - fio) (109) 

Theorem 6.6. Consider the instantaneously thermalized version of the Brownian energy 
process, with generator CJ^ . This process is dual to the process with generator C^^'^l^^ with 
duality function given by 

D{e,i)=J{D,{e„ii) (110) 
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D ie- e ) - ^^"^/^^ fill) 

Proof. By the duality result for the Brownian energy process, Theorem 16.21 we have for all 

therefore, 



L\';>D{e,0 = 4™^'''""i5(e,e) (112) 



, J- (m) , J- (m) ,dual 

lim(e*-^«j -id)D{e,^) = \im{e^^-j -id)D{e,0- (113) 

t—fOO t^oo 

The result then follows from the definition of the processes, together with lemma 16.31 and 
lemma 16.51 □ 



6.4 Limiting processes as m — > oo 

As it was done for the 2j-SEP, we study here the limiting behavior of the m-BEP process 
for large m. 

Theorem 6.7. Consider the process {z^.^^ '■ t >0} with generator L^™) and initial condition 
g and its dual {Ct"^^ : i > 0} with generator L^J^^i and initial condition G N*^. 
Suppose that, as oo, z^™) ^ z G and ^(""^ ^ ^ G N*^. Then: 

1. the process {z^J!^ ■ t > 0} converges to the process, {zt)t>o started from z, with gener- 
ator 

L = ^p{i,j)Lij 



- -2{z, - Zj) - ^) 



2. the process {'^^^^l converges to a system of independent random walkers {S,t)t>Q started 
from ^, with generator 

ijes 

Lfr' = mfic') - fio) + '^Cjifin - /(O) 

3. The two limiting processes {xt)t>o o,nd {£,t)t>o above are each other dual, with duality 
function 

Dix,0 = llxf 
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Proof. The proof of items 1. and 2. proceeds like in Theorem 14.31 For item 3. compare to 
example in section 13.51 □ 



6.5 Boundary driven process 

In this last section we consider the m-BEP process in contact at its boundary to energy 
reservoirs of the Ornstein-Uhlenbeck type. A duality result for the Brownian Momentum 
Process with reservoirs was already proven in [^. Here we generalize this result to the 
general Brownian energy process for arbitrary m G N. We start from the momentum process 
{{x(t)i^a) : i G 5, a = 1, . . . , m, t > 0} on a ladder graph with m levels and all levels at sites 
i S dS connected to a thermalizing Ornstein-Uhlenbeck process which parameter Tj, to be 
thought as the temperature. The generator reads 



m / 

i,jeS a,/9=l ^ 



d _ d 



™ ^2 Q 

+ y.y.'^iirr - ^^,o.J, (114) 

«ea<S a=l *>" ' 

If we now consider the induced process {zi{t) : i £ S,t >0} measuring the energy at each 
site via the map 

m 

E2 

a=l 

then, using the identities ([90|) . we find the generator 



i,jes 

+ E2r<fm|- + 2.,^)-2.<|- (115) 



ieds 



Introducing as usual a set d^S of sink sites and a bijection i ^ i^, which associate each 
boundary site i £ dS to a sink site ie S dgS, we have the following duality theorem: 

Theorem 6.8. Let {zt)t>o denote the boundary driven m-BEP with generator (jllSp . Then 
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{zt)t>o is dual to the process {(,t)t>o with generator 

LduaifiO = E ^'(^ J)2e.(2e, + m)(/(r'^) - fiO) + 2^,(2^. + m)ifien - fiO) 

+ E - /(C)) (116) 

ieds 

with duality function given by 

D(z, 0= T\ Tf^ TT 4' . ^ (117) 

Proof. The bulk part of the duahty function coincides with the one of Theorem 16. 2t the 
boundary part is easily checked with an explicit computation. □ 
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